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Abstract 

We present an approach to steady-state mesoscopic transport based on the 
maximum entropy principle formulation of nonequilibrium statistical mechan- 
ics. Our approach is not limited to the linear response regime. We show that 
this approach yields the quantization observed in the integer quantum Hall 
effect at large currents, which until now has been unexplained. We also pre- 
dict new behaviors of non-local resistances at large currents in the presence 
of dirty contacts. 
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In this Letter, we propose a non-perturbative general approach to nonlinear nonequi- 
librium steady-state transport in mesoscopic systems. Our work is based on the maximum 
entropy approach (MEA) to nonequilibrium statistical mechanics in which the den- 
sity matrix is found by by maximizing the information entropy of the system, subject to 
constraints which fix the expectation values of observables. Although the MEA should in 
principle be applicable to any nonequilibrium system, examples and explicit calculations 
have in practice been limited. In part this is because it is seldom possible to calculate 
the microstates of nonequilibrium systems. Moreover, it is difficult in general to determine 
whether the information entropy is equal to the thermodynamic entropy, as is needed, e.g., 
to identify the temperature in the density matrix derived from the MEA. The problem of 
steady-state mesoscopic transport, however, is uniquely well suited to this approach: the 
microstates can be calculated to obtain the exact density matrix, and recent work by Her- 
shfield @ allows us to identify the temperature in the density matrix. Here we show how 
the MEA can be used to calculate nonlinear current- voltage relations in mesoscopic devices. 

Our study was initially motivated by a very important but often neglected fact: the in- 
teger quantum Hall effect (IQHE) is exhibited even in systems driven by very large currents 
||. The IQHE can be viewed as a near-ideal manifestation of mesoscopic transport ||. An 
elegant explanation of the IQHE at low currents is given by the Landauer-Biittiker (LB) 
approach || to mesoscopic transport. However, this approach is fundamentally a linear 
response theory |||7[], and, as we show below, fails to yield the quantization observed at high 
currents ||. Hence the IQHE exhibits a 'simple' behavior (perfect quantization) far beyond 
the linear response regime. There exist other, general, approaches to nonequilibrium trans- 
port, such as various Green's function techniques PH12|. These lead to quite complicated 
calculations even in the linear regime, and it is not clear whether they can yield the IQHE 
in the nonlinear regime. Any comprehensive theory of nonlinear mesoscopic transport must 
be able to explain the extraordinary quantization of the IQHE at high currents. One of our 
most important results is to show that this can in fact be explained within the MEA. 

We consider a mesoscopic system consisting of a device (e.g. a Hall bar or a quantum 
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wire), to which M terminals, denoted by m (m = s, d, 1, 2, . . . , M — 2) are connected by 
straight leads long enough that evanescent modes emanating from the terminals decay to zero 
in the leads. A current / flows from the source s to the drain d. Electrons injected in states 
at one terminal can either be transmitted to another terminal, or reflected back. In either 
case, they lose their phase memory upon re-entering a terminal due to phase-randomizing 
scattering there. The electrons in the system are described by a complete orthogonal set 
of eigenstates \ip a ) with energies e a . Each eigenstate in general carries a net current i m>a 
from each terminal [ l3| . A particularly useful set of eigenstates for multi-terminal (M > 2) 



systems are the scattering states [0,|T4| iV'mnfc) usec ^ wnen the terminals are modeled as semi- 
infinite straight leads. The state \ipmnk) * s incoming into the device from terminal m; n and 
k denote the asymptotic wavenumber and subband index of the incoming wave. With a 



proper normalization ||14|| , this state's net current, i m ',mnk, at terminal m' is related to the 
incoming current i° mnk by i m >,mnk = ^mnk^m'm - J2n'k'tm'n'k',mnk), where t m ' n > k ^ mnk is the 
transition probability obtained from the scattering matrix in the iVwfc) representation. 

In the LB formalism || it is assumed that each terminal is held at a 'local chemical 
potential' /i m , so that electrons are injected into the device at each terminal with distribu- 
tions fmnk = l/[e^ em " fc ~ Mm * ) + 1]. For a two-terminal system at zero temperature and low 
voltages V = — /irf)/e, this gives a resistance R = h/(je 2 t). Here j is the number of 



occupied subbands and t is the total transmission probability at fj, s [fLq| . This resistance is 
quantized in the absence of backscattering (t = 1); this is the two-terminal version of the LB 
explanation of the low-current IQHE. However, when the voltage becomes greater than the 
subband spacing, the source injects electrons into the (j + l)st subband but the drain does 
not. Then, according to the LB approach, the resistance of an ideal system would lie between 
h/(je 2 ) and h/[(j + l)e 2 ]. The same conclusion is reached for the case of a multi-terminal 
system. (We point out that this argument has been invoked to explain the large-voltage 
failure of resistance quantization in quantum point contact experiments within the LB for- 
malism [jnj.) Yet in precision IQHE measurements, the voltage is many times greater than 
the subband spacing and the resistance is nonetheless highly quantized ||. Hartree inter- 



actions moves subbands together and cannot change this conclusion. Exchange-correlation 
interactions could restore quantization if they caused a large energy separation (of order eV) 
between occupied and unoccupied subbands. However, the exchange- correlation energies in 
the fractional quantum Hall effect, for example, are only of the order of 10~ 3 hu> c . Thus 
linear response theory is unable to explain the cleanest experiments in mesoscopics: the 
quantization found in IQHE systems at large currents. 

Let us now describe a general approach to nonlinear steady-state mesoscopic transport 
which resolves this. The thermodynamic variables which can be taken as known are the 
internal energy U and the particle number N. We add to these the net current I m at each 
terminal (so that I s = —Id = I and I m = 0, m 7^ s, d). Following the MEA U, we then 
maximize the information entropy Sj = —cJ2yP^ l n p 7 , subject to constraints on the average 
values of energy, particle number, and currents. Here c is an (as yet) unspecified constant, 
and p 1 is the probability that the system is in a microstate 7. This can be written as the 
matrix element p 1 = (7I/0I7) of the density matrix p. Averages of an operator A are given 
by (A) = Tr Ap /Trp. The constraints are imposed by requiring that (H) = U, (N) = N, 
and (J m ) = I m . Here H, and N are, respectively, the Hamiltonian and particle number 



operators, and I m is the net current operator in lead m [0 . The constrained maximization 
gives the density matrix 

p = exp [-/3(H - fiN - (1) 

rn 

In this expression p is the global chemical potential, associated with a global particle reser- 
voir, and the intensive variables £ m are Lagrangian multipliers associated with the constraints 
on the currents. Because of current conservation there are only M — 1 independent current 
constraints, so we may choose ^ = 0. The variable (3 is the product of c~ l and the vari- 
able conjugate to U. This density matrix has the general form which Hershfield recently 
showed exists quite generally in steady-state nonequilibrium systems; following his work 
we therefore identify (3 = 1/ksT, where kg is Boltzmann's constant and T the thermody- 
namic temperature. This identification also means that in this case the information and 



thermodynamic entropies are identical (with c = fee). 

This formal result can be more clearly written in terms of a complete set of single-particle 
eigenstates \ip a ) of H and I m . The above density operator then gives the following thermal 
occupancies of single-particle states: 

^ a 1 + exp [{3 (e Q - /I - £ m £ m i m ,a)] ' ^ 
For illustration, consider this result in terms of the scattering states of a two-terminal system. 
In this case we can drop the terminal index m, and understand that k > corresponds to 
states injected by the source and k < to states injected by the drain. In an ideal system 
with t n 'fc', ra fc — firin^k'k, these states carry currents i n k, and f n k = 1/ 
where £ = £ s . In the simplest case, with only one subband (n = 0) occupied, this is 
similar to the LB result; the combination fi + £iokt acts like an effective local chemical 
potential. However, with more than one subband occupied, or with nontrivial t's, or at large 
voltages, our current-constrained equilibrium occupations cannot be described in terms of 
local chemical potentials. In the general case states are occupied up to different energies in 
each subband. 

A voltmeter connected between the source and drain measures the work required to 
move a unit charge between them. In the LB formalism this voltage is simply (fi s — fid) / e. 
Let us find the corresponding result for our distribution. We use the representation given 
by the scattering states (V'mnfe)' which have occupancies f mn k- For clarity here we will 
present the results for a two-terminal device. The generalization to the multi-terminal case is 
straightforward and will be presented elsewhere. The steady-state condition and the absence 
of inelastic scattering within the device allow us to define thermodynamic potentials of the 
electron distribution, just as for an equilibrium system J20|| . For example, the equivalence 
in the present case between the information and thermodynamic entropies means that, as 
in equilibrium, here the Helmholtz free energy F = U — TS. The thermodynamic work 
5W done on the system at constant temperature is then equal to the change in free energy, 
5F = ji5N + J2m£,m SI m ', in the case of two terminals, this becomes SF = jj,5N + £ 51, where 
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I is the source-to-drain current, and £ = £ s . Varying either of the variables /x and £ generally 
changes the occupancy of states injected by both terminals, by Eq. (0). Let ry refer to either 
of the variables /x, £, and let 5I V be the change in net current when 77 is varied with the other 
variable held fixed. Similarly let SN^ be the change in the occupancy of scattering states 
injected by terminal m when r\ is varied (with N m = J2 n k fmnk the total particle number 
injected at terminal m), so that the change in free energy is SF V = fi(SN^ + SN]) + £ 8I n . 
We obtain the potentials V m at the terminals by interpreting this free energy change as the 
work done in adding electrons SN^ injected at each terminal against the voltage V m at the 
terminal. Thus, a change SN^ occurs at a cost in work of eSN^Vm. The total work is then 
e(SN^V s + 5N]Vd), and equating this to 5F V for r] = fx, £ gives two linearly independent 
equations 

E SNl(eV m -n)=^SP, (3) 

m=s,d 

which must be solved for the unknown terminal voltages V s and Vd- (In the M-terminal case, 
this becomes a set of M equations.) The resistance measured between source and drain is 
then R = (V s — Vd)/I. The 2x2 matrix 8N^ on the left-hand side of Eq. (^) is invertible, 
and the resulting potentials automatically are given relative to the global chemical potential 

The distribution f a in Eq. (0) has been written down in earlier work by Heinonen and 
Taylor [^TJ, who used it to study current distributions, and more recently by Ng [f[9[j . In 



these works it was argued that the lack of dissipation in a device without inelastic scattering 
permitted the ansatz of minimizing a free energy, subject to the current constraint. Here we 
have shown how this can be justified much more generally within the MEA, and the absence 
of dissipation makes it possible to determine all microstates. The second completely new 
point in the current work is our calculation of voltage from considerations of work. Ng, 
for example, simply assumed that the current-induced potential difference is proportional 
to the change in occupancies at the terminals. The validity of this assumption is not at all 
clear in, e.g., precision IQHE measurements where the Hall voltage is much greater than the 
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bulk Fermi energy. Furthermore, Ng failed to constrain particle number, and consequently 
predicted that even at small currents only states injected at the source (and not the drain) 
should be occupied. This appears unphysical. 

We will illustrate our approach with two examples, turning first to the resistance of 
an ideal two-terminal system. For simplicity, we we drop the terminal subscripts, and use 
eigenstates which satisfy periodic boundary conditions on a length L along the device. (This 
is only to choose a simple density of states; the final result does not depend on this particular 
boundary condition.) Suppose that current-carrying states have energies e n k = e n +h 2 k 2 /2m* 
and carry currents ehk/m*L. This can represent ID transport, or a parabolically confined 
Hall bar. The occupancies are, by Eq. @, f nk = f(e n k — £,i n k — t 1 ), where /(e) = l/(e /3e + 1). 
Then f n k is symmetric about k = £ = £e/hL, and we define e n k = e n + h 2 (k — £) 2 /2m* 
so that e n k — fi — £,i n k = ?nk — U, where Jx = li + h 2 ^ 2 /2m*. The electron number is 
N = Z)n(2vr/L) J dk f(e n k — Ll), and the total current is / = J^ni^/L) J dki n kf(e n k — fi)- We 
convert the integrals over k to integrals over energy e nk , and obtain iV = J2 n IT dep n (e)f(e — 
jl) where p„(e) = (L / n){2h 2 (£ — e n )/m*]~ 1//2 is the ID density of states. Similarly, / = 
e%/{m*L)^ X]n(27r/-^) / dk f{e n k — fi>), so I = ehN^/m* L. It is then simple to calculate 8I V 
and 5N^, in Eq. (^), with the integrals for 8N^ over k restricted to k > (k < 0) for 
m = s(d). The resulting expression for the voltage difference V s — Va obtained by inverting 
Eq. @ is simplified by the symmetry of f n k about £, and we find 

« = «-^ = ?(??e^t)" 1 - « 

If /i exceeds only the band minima of the first j subbands (or Landau levels), then at 
zero temperature R = h/je 2 . Finite-temperature corrections are exponentially small. This 
exact result is true regardless of the size of the voltage or current. This is perhaps our 
most important result; ours is the first mesoscopic transport theory which can explain the 
extremely accurate quantization seen in the IQHE far from the linear response regime. (We 
neglect the breakdown which occurs in the IQHE at very large current densities when other 



dissipative mechanisms turn on p2| .) If, for I > 0, there are no states with k < occupied, 



R is not quantized; this appears to be the case with quantum point contacts at large currents 
fT6| . We have numerically studied non-parabolic energies e n k and multi-terminal systems 
and find in these cases the accuracy of the quantization is limited only by the numerical 
accuracy, so long as there are states with k < occupied. 

As a second example, we consider a system with a 'dirty' source (a source with backscat- 
tering). The LB formalism involves only Fermi surface properties H0. As a consequence 
23fl , even in the presence of such 'dirty contacts', in the LB approach all Hall conductances 



are quantized and all longitudinal conductances are zero — provided that no two dirty con- 
tacts are adjacent to one another. This in fact can only be true in the linear response 
regime. At finite current J, the net current at each terminal current involves an integral 
over transmission probabilities, and resistances need then not be quantized. Consider as a 
simple example the four-terminal resistances in a system with a dirty source in the presence 
of a magnetic field, with only the lowest subband occupied. We find that only at very small 
currents do the resistances attain their ideal values (this is illustrated in Fig. [I]). Even the 
LB approach, if applied naively beyond linear response, gives deviations in some resistances. 
To inject a given current despite back-reflection at the source requires /i s to be greater than 
its value in the ideal case. Then resistances between the source and the other terminals 
differ from the ideal values. Resistances not involving the source are still ideal (quantized 
or zero) in the LB formalism. In our approach, even the latter are non-ideal at finite cur- 
rents (see Fig. [I]) because the occupancies of electrons injected at terminal m depend on the 
transmission from all other terminals into m [Eq. (0)]. 

In the MEA observables enter the formalism as constraints. Consequently, we have 
included the presence of a current /, driven by a current source, as a constraint on the net 
current. The ability of this to describe the IQHE at large currents is not trivial and argues, 
we believe, for its validity. We note that the LB formalism can also be obtained from the 
MEA if the current source is assumed to constrain the particle number N m injected at each 
terminal, rather than the current. These constraints are imposed by Lagrangian multipliers 
fi m , and the occupancies which result are precisely the LB distributions f^nk- At l° w 
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currents the use of local chemical potentials can be justified using linear response theory, 
viewing the potential difference (or the associated electric field) as driving the current. 
This cannot be extended to high currents. Notice that in the MEA the driving force need 
not be represented by an operator in the Hamiltonian. Instead, the result of the driving 
(here, the current) enters as a constraint. In the MEA the LB distribution would arise if 
a current source could be thought of as an entity that controls particle number instead of 
current. One might suppose that this models a voltage source instead of a current source. 
If so, then the I-V curve at large currents and voltages would depend on whether voltage 
or (as is usual) current is applied ||19|| . (In the linear regime, both approaches give the 
same result.) This is possible in principle, since voltage differences correspond to work, 
which is not a thermodynamic state function. As we now show, this appears not to be the 
case, and distributions of the form Eq. (0) should in fact be expected in a steady-state 
dissipationless system. Consider an arbitrarily long ideal device in which electrons flow in 
their steady-state distributions with a net current. In the device the Hamiltonian (including 
electron-electron-electron interactions and electron-phonon interactions) is translationally 
invariant (neglecting Umklapp processes) and preserves the distributions. In the frame of 
reference moving with velocity v such that the net current vanishes, the electrons are then in 
equilibrium at some chemical potential //. Hence the occupancy of a state with energy e' nk is 
in this frame given by f(e' n k — //). By a Galilean transformation, e' nk — y! = e n k + hkv — fi = 
tnk + (.ink — where proportionality constants have been absorbed in £, and \i differs from 
// by a constant. Since the occupancy of a state is the same in each frame, the distributions 
in the stationary frame are thus f{e n k — f> ~ £,ink)- 

The distributions Eq. (0) lead to other observable phenomena. For example, they lead to 
dissipationless deviations in quantization in the IQHE when states in different subbands are 



mixed by short-range elastic scatterers [24]. This might explain recent observations in high- 



quality Si samples |{25|. We have also used our formalism to explain the I — V characteristics 



of quantum point contacts and the lack of current saturation at high voltages |3B| . 
The approach we have presented includes nonlinear effects due to the current-dependent 
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electron distributions. At higher currents, other nonlinearities arise from distortions of the 
electron wavefunctions by the resulting electric field. This field is due to electron-electron 
interactions, which can easily be included in our approach at the Hartree level. In preliminary 
numerical calculations this causes no qualitative change in the picture. 

This work was supported in part by the UCF Division of Sponsored Research. 
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FIGURES 

FIG. 1. The resistance R s d,id measured between terminal 1 and the drain when current flows 
from source to drain in a four-terminal IQHE system with a dirty source. The device is a cross, 
with terminals 1 and 2 on opposite sides of the current flow. In an ideal system, this resistance is 
h/e 2 at zero temperature. The curves are at temperatures T = 0.0025, 0.005, 0.01, 0.05, and 0.1 
in units of huj c /kB- The zero-current Fermi energy is 1.32 hcj c . Insert: The reflection probability 
r at the source for the lowest subband. 
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